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ABSTRACT 

The structure of three-jet e"'"e~ — >• bbg events has been studied using hadronic Z° decays 
recorded in the SLD experiment at SLAC. Three-jet final states were selected and the 
CCD-based vertex detector was used to identify two of the jets as b or 6; the remaining 
jet in each event was tagged as the gluon jet. Distributions of the gluon energy and 
polar angle with respect to the electron beam were measured over the full kinematic 
range, and used to test the predictions of perturbative QCD. We find that beyond- 
leading-order QCD calculations are needed to reproduce the features seen in the data. 
The energy distribution is sensitive to an anomalous b chromomagnetic moment k, at 
the bbg vertex. We measured k to be consistent with zero and set 95% C.L. limits on 
its value, -0.06 < k < 0.04. 
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1. Introduction 



The observation of e^e~ annihilation into final states containing three hadronic jets, 
and their interpretation in terms of the process e^e" —>■ qqg provided the first 
direct evidence for the existence of the gluon, the gauge boson of the theory of strong 
interactions, Quantum Chromodynamics (QCD). In subsequent studies the jets were 
usually energy ordered, and the lowest-energy jet was assigned as the gluon; this is 
correct roughly 80% of the time, but preferentially selects low-energy gluons. If the 
gluon jet can be tagged explicitly, event-by-event, the full kinematic range of gluon 
energies can be explored, and more detailed tests of QCD can be performed 0. Due 
to advances in vertex-detection this is now possible using e"'"e~ bbg events. The large 
mass and relatively long lifetime, ~ 1.5 ps, of the leading B hadron in 6-quark jets 
lead to decay signatures that distinguish them from lighter-quark [u, d, s or c) and 
gluon jets. We used our charge-coupled-device (CCD)-based vertex detector (VXD) ^ 
to identify in each event the two jets that contain the B hadrons, and hence to tag the 
gluon jet. This allowed us to measure the gluon energy and polar-angle distributions 
over the full kinematic range. 

Additional motivation to study the bbg system has been provided by measure- 
ments involving inclusive bb decays. Several early determinations of -R^ = 
T{Z^ -^bb)/T{Z^ ^qq) differed from Standard Model (SM) expectations at the few 
standard deviation level. More recently it has been noted that the LEP measurement 
of the 6-quark forward-backward asymmetry, A^p^, lies roughly 2.5 standard devia- 
tions below the SM expectation. Since one expects new high-mass-scale dynamics 
to couple to the massive third-generation fermions, these measurements have aroused 
considerable interest and speculation. We have therefore investigated in detail the 
strong-interaction dynamics of the 6-quark. We have compared the strong coupling 
of the gluon to 6-quarks with that to light- and charm-quarks 0, as well as tested 
parity (P) and charge©parity (CP) conservation at the bbg vertex 0]. We have also 
studied the structure of bbg events, via the distributions of the gluon energy and polar 
angle with respect to (w.r.t.) the beamline [§], using the JADE algorithm [^j for jet 
definition. 

Here we update the bbg structure measurements using a data sample more than 3 
times larger than in our earher study, and recorded in the upgraded vertex detector. 
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which allowed us to improve significantly the gluon-jet tagging. In addition we extended 
our study to include the Durham, Geneva, E, EO and P algorithms |ffU[ to define jets. 



and compared these results with perturbative QCD predictions. This constitutes a 
more detailed test of QCD and enabled us to study systematic effects arising from the 
jet definition algorithm. 

Furthermore, we have used these data to study possible deviations from QCD in 
the form of radiative corrections induced by new physics, in terms of the 6-quark 
chromomagnetic moment. In QCD this is induced at the one-loop level and is of order 



as/vr. A more general bbg Lagrangian term with a modified coupling IT^] may be 
written: 

C''^ = g.bTaii, + '-^{k - m75)}&GJ^, (1) 

where k and k parametrize the anomalous chromomagnetic and chromoelectric mo- 
ments, respectively, which might arise from physics beyond the SM. The effects of the 
chromoelectric moment are sub-leading w.r.t. those of the chromomagentic moment, 
so for convenience we set k to zero. A non-zero k would be observable as a modifica- 



tion of the gluon energy distribution in hhg events relative to the standard QCD 
case. We have used our precise measurements of the gluon energy distributions to set 
the most stringent limits on n. 

2. bbg Event Selection 

We used hadronic decays of bosons produced by e+e^ annihilations at the SLAC 
Linear Collider (SLC) and recorded in the SLC Large Detector (SLD) [T^. The criteria 
for selecting hadronic Z*^ decays and the charged tracks used for fiavor-tagging are 



described in 0, 0- Three-jet events were selected using iterative clustering algorithms 
applied to the set of charged tracks in each event. We used in turn the JADE, Durham, 
E, EO, P and Geneva algorithms. The respective scaled- invariant-mass, Ucut-, values were 
chosen to maximise the statistical power of the measurement, while keeping systematic 
errors small; they are shown in Table |l|. 

Events classified as 3-jet states were retained if all three jets were well contained 
within the barrel tracking system, with polar angle |cos^jet| ^ 0.80. In addition, 
in order to select planar 3-jet events, the sum of the angles between the jet axes 
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Jet algorithm 


Ucut value 


7^ 3-jet events 


efficiency 


T A r>p 


u.uzo 


K7'^A 1 
( 04i 


1 9 9% 
IZ.Z/o 


Durham 


0.0095 


46432 


12.1% 


E 


0.0275 


66848 


11.7% 


EO 


0.0275 


54163 


11.2% 


P 


0.02 


60387 


12.0% 


Geneva 


0.05 


40895 


12.8% 



Table 1: Number of selected 3-jet events, and gluon-jet tagging efficiency (see text), 
for each algorithm. The statistical error on the efficiency is roughly 0.04%. 



was required to be between 358 and 360 degrees. From our 1996-98 data samples, 
comprising roughly 400,000 hadronic decays, the numbers of selected events are 
shown in Table |l[ In order to improve the energy resolution the jet energies were 
rescaled kinematically according to the angles between the jet axes, assuming energy 
and momentum conservation and massless kinematics. The jets were then labelled in 
order of energy such that Ei > E2 > E^. 

Charged tracks with high quality information in the VXD as defined in Q were 
used to tag hbg events. For each such track we examined the impact-parameter, d, 
w.r.t. the interaction point (IP). The resolution on d is given by ad =7.7©29/psin^''^ 9 
/im in the plane transverse to the beamline, and 9.6©29/psin^''^ 6' /im in any plane 
containing the beamline, where p is the track momentum in GeV/c, and 9 the polar 
angle, w.r.t. the beamline. 



Jets containing heavy hadrons were tagged using a topological algorithm ap- 
plied to the set of tracks associated with each jet. A track density function was calcu- 
lated, and a region of high total track density well separated from the IP was identified 
as a vertex from the decay of a heavy hadron. For each vertex, the pt-corrected invari- 



ant mass |T^, Mpj, was calculated using the set of tracks associated with the vertex, 
assuming the charged pion mass, and the vertex axis, defined to be the vector from 
the IP to the reconstructed vertex position. Fig. |l] shows the Mp^ distributions sepa- 
rately for vertices found in jets 1, 2 and 3 using, for illustration, the JADE algorithm; 
results using the other algorithms (not shown) are qualitatively similar. The simu- 
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lated contributions from true b, c, light and gluon jets are indicated; the c, hght and 
gluon jets populate predominantly the region Mp^ < 2 GeV/c^. Events were retained 
in which exactly two jets contained such a vertex, and at least one of them had Mp^ 
> 2 GeV/c^. In order to suppress events in which a single 5-hadron decay gave rise 
to two reconstructed vertices, the cosine of the angle between the two vertex axes was 
required to be less than 0.7, and the distance between the vertices, projected in a plane 
perpendicular to the beamline, was required to exceed 0.12 cm. Roughly 1.1% of the 
event sample was rejected by these cuts. In each selected event the jet without a vertex 
was tagged as the gluon jet. 

For each algorithm, the number of tagged jets is shown in Table 0; also shown, 
in Table is the efficiency for tagging the gluon jet correctly in true bbg events, 
which was calculated using a simulated event sample generated with JETSET 7.4 |]15 



with parameter values tuned to hadronic e^e~ annihilation data |T^, combined with a 
simulation of 5-decays tuned to T(4S) data and a simulation of the detector. For 
the JADE algorithm, for example, the efficiency peaks at about 15% for 18 GeV gluons. 
Below 18 GeV the efficiency falls, to as low as 3%, since lower-energy gluon jets are 
sometimes merged with the parent &-jet by the jet-finder. Above 18 GeV the efficiency 
falls, to as low as 5%, since at higher gluon energies the correspondingly lower-energy 
6-jets are more difficult to tag, and there is also a higher probability of losing a jet 
outside the detector acceptance. Results for the other algorithms are qualitatively 
similar. The systematic error associated with the tagging efficiency was small and was 
explicitly taken into account by the procedure for estimating systematic uncertainties 
that is described in Section 3. 

For each algorithm the tagging purities, defined as the fraction of selected 3-jet 
events in which the two vertices were found in the two jets containing the true B 
hadrons, are listed by gluon-jet number in Table We formed the distributions of two 
gluon-jet observables, the scaled energy Xg = 2Egi^on/ Vs, and the polar angle w.r.t. 
the beamline, 6g. For illustration, for the JADE algorithm the distributions are shown 
in Fig. ^ the simulation is also shown; it reproduces the data. Results for the other 
algorithms (not shown) are qualitatively similar. 

The backgrounds were estimated using the simulation and are of three types: non- 
66 events; bb but non-bbg events; and mis-tagged true bbg events. Their contributions 
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Figure 1: The Mpt distributions for vertices found in selected 3-jet events, defined using 
the JADE algorithm, labelled according to jet energy (dots); errors are statistical. 
Histograms: simulated distributions for different jet flavors. Events were selected by 
requiring that at least two jets contain a vertex, at least one of which must satisfy Mpt 
> 2 GeV/c2 (see text). 
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Figure 2: Raw measured distributions of (a) Xg and (b) cos9g (dots) defined using the 
JADE algorithm; errors are statistical. Histograms: simulated distributions including 
background contributions. 
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JADE 


Durham 


E 


jet label 


# jets 


purity (%) 


# jets 


purity (%) 


# jets 


purity (%) 


3 


4349 


98.0 


2952 


97.0 


5246 


97.5 


2 


740 


90.2 


890 


92.4 


1007 


85.4 


1 


150 


71.0 


138 


73.4 


148 


70.8 




EO 


P 


Geneva 


jet label 


# jets 


purity (%) 


# jets 


purity (%) 


# jets 


purity (%) 


3 


4027 


98.0 


4654 


98.0 


3491 


93.9 


2 


692 


90.2 


795 


90.7 


692 


86.7 


1 


151 


70.7 


155 


72.1 


181 


63.4 



Table 2: Tagging purities (see text). 

are shown in Fig. § for the JADE case. Results for the other algorithms (not shown) 
are qualitatively similar. For each algorithm, the non-66 events make up roughly 1% 
of the selected sample and are dominated by ccg events; roughly 70% of these had the 
gluon jet correctly tagged, and the remainder comprises events in which the gluon split 
into a cc or 66, which yielded a real secondary vertex in the 'wrong' jet. Mis-tagged 
events, in which the gluon jet was mis-tagged as a 6 or 6-jet and one of the true 6- or 
6-jets enters into the measured gluon distributions, comprise roughly 3% of the sample; 
roughly two thirds of these events contain a gluon splitting into cc or 66. These two 
backgrounds are negligible except in the highest Xg bin. 

For all algorithms the dominant background {eg. for JADE, roughly 16% of the 
sample) is formed by 66 but non-bbg events. These are true 66 events that were not 
classified as 3-jet events at the parton level, but were reconstructed and tagged as 3-jet 
bbg events in the detector using the same jet algorithm and ycut value. In a parton-level 
2-jet event this can arise from the broadening of the particle flow around the original 
6 and 6 directions due to hadronization and weak decay; in particular, the relatively 
high-transverse-momentum 5-decay products can cause the jet-finder to reconstruct 
a 'fake' third jet, which is almost always assigned as a (low-energy) gluon jet. In 
addition, an event classified as 4-jet at the parton level may, due to the overlap of their 
hadronization products, have two of its jets combined in the detector by the jet-finder. 
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In this case the combined jet is usually tagged as a gluon jet . Since the calculations 
with which we compare below are not reliable for 4-jet events, we consider such events 
to be a background. 



3. Correction of the Data 

For each algorithm, the distributions were corrected to obtain the true gluon distri- 
butions by applying a bin-by-bin procedure: = C(X) (^'^'^'"(X) - 
B{X)), where X = or cosOg, Z}'™"'(X) is the raw distribution, B{X) is the back- 
ground contribution, and C{X) = Dj7c("^)/-^mc"("^) ^ correction that accounts 
for the efficiency for accepting true bbg events into the tagged sample, as well as for 
bin-to-bin migrations caused by hadronization, the resolution of the detector, and bias 
of the jet-tagging technique. Here D*^^q{X) is the true distribution for MC-generated 
bbg events, and D]y^™"'(X) is the resulting distribution after full simulation of the de- 
tector and application of the same analysis procedure as applied to the data. 

The fully-corrected distributions are shown in Figs. ^, |, ^, ^ |^, and ^. Since, in an 
earlier study 0, we verified that the overall rate of bbg-event production is consistent 
with QCD expectations, we normalised the gluon distributions to unit area and we 
study further the distribution shapes. In each case the peak in Xg is a kinematic 
artefact of the jet-finding algorithm, which ensures that gluon jets are reconstructed 
with a non-zero energy, and it depends on the ycut value. The cos^^ distributions are 
very nearly fiat, in contrast to the 1 + cos^ 6 behavior for quark jets. 

We have considered sources of systematic uncertainty that potentially affect our 
results. These may be divided into uncertainties in modelling the detector and uncer- 
tainties in the underlying physics modelling. To estimate the first case we systemat- 
ically varied the track and event selection requirements, as well as the track-finding 
efficiency [^, the momentum and dip angle resolution, and the probability of finding 



a fake vertex in a jet. In the second case parameters used in our simulation, relating 
mainly to the production and decay of charm and bottom hadrons, were varied within 
their measurement errors HTSf. For each variation the data were recorrected to de- 



rive new Xg and cosOg distributions, and the deviation w.r.t. the standard case was 
assigned as a systematic uncertainty. Although many of these variations affect the 
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Figure 3: Corrected distributions of (a) Xg and (b) cos6'g (dots) defined using tlie 
JADE algoritfim; tlie error bars represent tlie sum in quadrature of tlie statistical and 
systematic errors. Perturbative QCD predictions (see text) are shown as lines joining 
entries plotted at the respective bin centers. 
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Figure 4: Corrected distributions of (a) Xg and (b) cos9g (dots); defined using the 
Durham algorithm; the error bars represent the sum in quadrature of the statistical 
and systematic errors. Perturbative QCD predictions (see text) are shown as lines 
joining entries plotted at the respective bin centers. 
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Figure 5: Corrected distributions of (a) Xg and (b) cos6'g (dots); defined using the 
E algoritlim; the error bars represent the sum in quadrature of the statistical and 
systematic errors. Perturbative QCD predictions (see text) are shown as hues joining 
entries plotted at the respective bin centers. 
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Figure 6: Corrected distributions of (a) Xg and (b) cos6'g (dots); defined using the 
EO algorithm; the error bars represent the sum in quadrature of the statistical and 
systematic errors, erturbative QCD predictions (see text) are shown as lines joining 
entries plotted at the respective bin centers. 
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Figure 7: Corrected distributions of (a) Xg and (b) cos6'g (dots); defined using the 
P algoritlim; the error bars represent the sum in quadrature of the statistical and 
systematic errors. Perturbative QCD predictions (see text) are shown as hues joining 
entries plotted at the respective bin centers. 
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Figure 8: Corrected distributions of (a) Xg and (b) cos9g (dots); defined using tlie 
Geneva algorithm; the error bars represent the sum in quadrature of the statistical and 
systematic errors. Perturbative QCD predictions (see text) are shown as lines joining 
entries plotted at the respective bin centers. 
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overall tagging efficiency, most had little effect on the energy or polar angle depen- 
dence, and no variation affects the conclusions below. The largest contributions to the 
error arose from the measurement uncertainties on the probability for gluon splitting 
into bb (which dominates around Xg ~ 0.5) or cc (which dominates for Xg > 0.7). 

All uncertainties were conservatively assumed to be uncorrelated and were added 
in quadrature in each bin of Xg and cos6g. In any bin the systematic error is typically 
much smaller than the statistical error. The data points with their total error bars are 
shown in Figures ^ |, |^, |^, ||; the data are listed in Tables ^ and ^. In addition, 
as cross-checks, for each algorithm the i/cut value and the Mp^ cut were varied around 
their respective default values; in no case did our conclusions change. 

4. Comparison with QCD Predictions 

We compared the data with perturbative QCD predictions for the respective jet algo- 
rithm and Ucut value. We calculated leading-order (LO) and next-to-LO (NLO) pre- 
dictions using JETSET. We also derived these distributions using the 'parton shower' 
(PS) implemented in JETSET; this is operationally equivalent to a calculation in which 
leading and next-to-leading Ini/c terms are partially resummed to all orders in a^. In 
physical terms this allows events to be generated with multiple orders of parton radi- 
ation, in contrast to the maximum number of 3 (4) partons allowed in the LO (NLO) 
calculations, respectively. Configurations with > 3 partons are relevant to the ob- 
servables considered here since they may be resolved as 3-jet events by the jet-finding 
algorithm. These predictions are shown in Figs. ^ ||, ^ |, |^ and ||. 

In the case of the cos^g distributions the three calculations are indistinguishable 
and they reproduce the data. For clarity we show in Figs. ||, ||, H, 0, and | only the 
PS calculations. We conclude that the cos6'g observable is insensitive to the details of 
higher order soft parton emission. 

In the case of Xg, for the JADE, E, EO and P algorithms the LO calculation repro- 
duces the main features of the shape of the distribution, but it yields too few events in 
the region 0.2 < Xg < 0.5, and too many events for < 0.1 and Xg > 0.6. The NLO 
calculation shows qualitatively similar behavior, although it reproduces the data no- 
ticeably better, especially for Xg > 0.6. In the case of the JADE, E and EO algorithms 
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1/N dN/dXg 


Xg range 




1 ) 1 1 V 1 1 >n 1 1 


F 


0.0- 


-0.1 


697±0 055±0 001 


000±0 000±0 000 


712±0 045±0 001 

\J m 1 J. ^ 1 \J m\J^I.\J 1 \J m\J\J ^ 


0.1- 


-0.2 


2.461±0.093±0.002 


1.093±0.073±0.005 


2.634±0.091±0.005 


0.2- 


-0.3 


2.021±0.074±0.003 


2.209±0.106±0.011 


2.111±0.075±0.010 


0.3- 


-0.4 


1.544±0.064±0.006 


1.982±0.093±0.018 


1.574±0.065±0.015 


0.4- 


-0.5 


1.158±0.056±0.007 


1.606±0.083±0.021 


1.111±0.055±0.019 


0.5- 


-0.6 


0.783±0.047±0.008 


1.196±0.072±0.015 


0.684±0.045±0.024 


0.6- 


-0.7 


0.559±0.042±0.008 


0.814±0.060±0.012 


0.493±0.041±0.029 


0.7- 


-0.8 


0.371±0.038±0.009 


0.501±0.049±0.010 


0.337±0.039±0.033 


0.8- 


-0.9 


0.287±0.038±0.012 


0.385±0.048±0.011 


0.231±0.038±0.032 


0.9- 


-1.0 


0.119±0.028±0.013 


0.214±0.043±0.013 


0.114±0.026±0.020 


Xg range 




p 


( --on f^^rsi 
vjrcii" vcL 


0.0- 


-0.1 


620-1-0 050-1-0 001 


1 362-1-0 092-1-0 001 


912-1-0 096-1-0 004 


0.1- 


-0.2 


2.369±0.092±0.002 


2.627±0.092±0.003 


2.582±0.130±0.008 


0.2- 


-0.3 


2.007±0.076±0.004 


1.779±0.064±0.004 


1.768±0.085±0.007 


0.3- 


-0.4 


1.509±0.065±0.006 


1.317±0.054±0.006 


1.304±0.070±0.007 


0.4- 


-0.5 


1.176±0.058±0.007 


1.033±0.049±0.007 


1.049±0.061±0.006 


0.5- 


-0.6 


0.817±0.050±0.007 


0.659±0.041±0.007 


0.713±0.052±0.007 


0.6- 


-0.7 


0.604±0.045±0.008 


0.511±0.039±0.008 


0.596±0.049±0.005 


0.7- 


-0.8 


0.411±0.041±0.009 


0.367±0.036±0.009 


0.433±0.047±0.008 


0.8- 


-0.9 


0.351±0.044±0.014 


0.268±0.036±0.012 


0.386±0.055±0.015 


0.9- 


-1.0 


0.137±0.033±0.016 


0.076±0.016±0.008 


0.259±0.064±0.028 



Table 3: Fully-corrected differential cross-sections for hard gluon-jet production as a 
function of jet energy Xg. The first error is statistical, the second systematic 
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1/N dN/dcos^g 


C()S/9y 


If Hi 1^0 


JADE 


1 ) 1 1 V 1 1 ^^ 1 1 

J ' I LJ. IJLfTllJLl 


E 


-0.80- 


-0.64 


503±0 035±0 004 


471±0 043±0 006 


497±0 036±0 010 


-0.64- 


-0.48 


0.512±0.027±0.003 


0.471±0.032±0.005 


0.533±0.028±0.009 


-0.48- 


-0.32 


0.484±0.026±0.003 


0.504±0.032±0.006 


0.470±0.025±0.009 


-0.32- 


-0.16 


0.482±0.027±0.003 


0.483±0.033±0.006 


0.467±0.025±0.009 


-0.16- 


-0.0 


0.530±0.028±0.003 


0.515±0.033±0.006 


0.528±0.026±0.009 


0.0- 


-0.16 


0.471±0.027±0.003 


0.481±0.032±0.005 


0.476±0.025±0.008 


0.16- 


-0.32 


0.494±0.027±0.003 


0.481±0.032±0.005 


0.472±0.025±0.008 


0.32- 


-0.48 


0.537±0.028±0.003 


0.509±0.033±0.005 


0.535±0.027±0.008 


0.48- 


-0.64 


0.520±0.027±0.003 


0.604±0.035±0.006 


0.539±0.028±0.009 


0.64- 


-0.80 


0.466±0.035±0.003 


0.480±0.044±0.005 


0.483±0.036±0.009 


K^KJoU g 




F,0 


p 


vjrciic vd 


-0 SO- 


-0 64 


4Q6-I-0 035-1-0 003 


536-1-0 035-1-0 004 


496-1-0 039-1-0 003 


-0.64- 


-0.48 


0.508±0.028±0.003 


0.511±0.027±0.003 


0.490±0.033±0.004 


-0.48- 


-0.32 


0.497±0.027±0.003 


0.486±0.025±0.004 


0.511±0.034±0.004 


-0.32- 


-0.16 


0.477±0.027±0.003 


0.463±0.025±0.003 


0.507±0.035±0.004 


-0.16- 


-0.0 


0.510±0.028±0.003 


0.508±0.026±0.003 


0.492±0.036±0.003 


0.0- 


-0.16 


0.469±0.027±0.003 


0.477±0.025±0.003 


0.547±0.037±0.004 


0.16- 


-0.32 


0.497±0.027±0.003 


0.490±0.026±0.003 


0.477±0.034±0.004 


0.32- 


-0.48 


0.519±0.028±0.003 


0.525±0.027±0.003 


0.501±0.035±0.004 


0.48- 


-0.64 


0.533±0.028±0.003 


0.538±0.027±0.003 


0.550±0.034±0.004 


0.64- 


-0.80 


0.493±0.037±0.004 


0.466±0.034±0.004 


0.429±0.038±0.003 



Table 4: Fully corrected differential cross-sections for hard gluon-jet production as a 
function of jet polar angle, cos9g. The first error is statistical, the second systematic 
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the PS calculation provides the best description of the data across the full Xg range, 
although it tends to underestimate the height of the kinematic peak; in the case of 
the P algorithm the PS calculation is slightly worse than the NLO calculation. These 
results suggest that the data are sensitive to multiple orders of parton radiation, the 
details of which need to be included in the perturbative QCD calculation. This is in 
agreement with our earlier inclusive measurement of jet energy distributions (for the 
JADE algorithm only) using flavor- inclusive Z° decays 

In the case of Xg defined using the Geneva algorithm (Fig. ^), there are clear dif- 
ferences among the three calculations, but the NLO calculation reproduces the data 
best. Finally, in the case of the Durham algorithm (Fig. the differences among the 
three calculations are relatively small, and both the NLO and PS calculations provide 
a good description of the data. This is consistent with the original motivation for the 
Durham algorithm [^], which was explicitly designed to yield a jet structure that is 
relatively insensitive to the presence of additional soft partons. 

We conclude that perturbative QCD in the PS (JADE, Durham, E, EO algorithms) 
or NLO (P, Durham, Geneva algorithms) approximation reproduces the gluon distribu- 
tions in bbg events. However, it is interesting to consider the extent to which anomalous 
chromomagnetic contributions are allowed by the data. The Lagrangian represented 
by Eq. |l] yields a model that is non-renormalizable. Nevertheless tree-level predictions 
can be derived [|1T| and used for a 'straw man' comparison with QCD. For each jet algo- 
rithm, in each bin of the Xg distribution, we parametrised the leading-order effect of an 
anomalous chromomagnetic moment and added it to the PS calculation to arrive at an 
effective QCD prediction including the anomalous moment at leading-order. A min- 
imization fit was performed to the data with k as a free parameter. The corresponding 
K and values are shown in Table |^. In all cases n is consistent with zero. For each 
algorithm the confidence level of the fit is smaller than the confidence level based on 
the foi' the comparison with the standard PS calculation. We conclude that our 
data show no evidence for any beyond-QCD effects, and we set 95% confidence-level 
(C.L.) limits on k; these are shown in Table |^. Since the results are highly correlated, 
we quote best limits on k using the JADE algorithm, yielding —0.058 < k, < 0.043 at 
the 95% C.L. 
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Jet algorithm 


K 




yo/o o.ij. limits 


T A mn" 


— u.uuo m u.uzu 


1 Q 
io.y 


— U.UOo K <^ U.U^tO 


Durham 


0.020 ± 0.043 


21.8 


-0.065 <K< 0.106 


E 


-0.005 ± 0.028 


13.6 


-0.060 <K< 0.050 


EO 


-0.006 ±0.027 


15.6 


-0.060 < K< 0.047 


P 


-0.002 ±0.025 


31.7 


-0.052 <K< 0.047 


Geneva 


-0.006 ±0.026 


24.4 


-0.056 <K< 0.045 



Table 5: Best-fit k, values and 95% C.L. limits. 

5. Conclusion 

In conclusion, we used the precise SLD tracking system to tag the gluon in 3-jet 
e~^e^ ^ Z° — s> bbg events. We studied the structure of these events in terms of 
the scaled gluon energy and polar angle, measured across the full kinematic range. 
We compared our data with perturbative QCD predictions and found that beyond-LO 
QCD contributions are needed to describe the energy distribution. We also investigated 
an anomalous 6-quark chromomagnetic moment, k, which would affect the shape of the 
energy distribution. We set 95% c.l. limits of —0.06 < k, < 0.04. These results are 
consistent with, more precise than, and supersede those in our earlier publication [Q. 

We thank the personnel of the SLAG accelerator department and the technical 
staffs of our collaborating institutions for their outstanding efforts on our behalf. We 
thank A. Brandenburg, P. Uwer and T. Rizzo for many helpful discussions. 
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